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ABSTRACT 

We examine the case where a circumstellar medium around a supernova is 
sufficiently opaque that a radiation dominated shock propagates in the circum- 
stellar region. The initial propagation of the shock front into the circumstellar 
region can be approximated by a self-similar solution that determines the ra- 
diative energy in a shocked shell; the eventual escape of this energy gives the 
maximum luminosity of the supernova. If the circumstellar density is described 
by p = Dr~^ out to a radius R^, where D is a constant, the properties of the 
shock breakout radiation depend on and Rd = nDvsh/c, where k is the opac- 
ity and Vsh is the shock velocity. If i?^ > Rd, the rise to maximum light begins at 
~ Rd/vsh', the duration of the rise is also ~ Rd/vgh] the outer parts of the opaque 
medium are extended and at low velocity at the time of peak luminosity; and a 
dense shell forms whose continued interaction with the dense mass loss gives a 
characteristic flatter portion of the declining light curve. If R^ < Rd-, the rise 
to maximum light begins at Rw/vgh', the duration of the rise is R^/vghRd, the 
outer parts of the opaque medium are not extended and are accelerated to high 
velocity by radiation pressure at the time of maximum luminosity; and a dense 
shell forms but does not affect the light curve near maximum. We argue that 
SN 2006gy is an example of the first kind of event, while SN 2010gx and related 
supernovae are examples of the second. 

Subject headings: circumstellar matter — shock waves — supernovae: general — 
supernovae: individual (SN 2006gy) 



INTRODUCTION 



Supernova shock breakouts from normal massive stars typically give rise to X-ray /ultraviolet 
bursts with a timescale < 10^ s ( Klein fc Chevalier|l978 Falk|l978 Ensman fc Burrows|l992 
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Matzner &: McKee|1999 Nakar fc Sari|20l"0 ). The timescale is generally determined by R^jc^ 
where i?^, is the stellar radius because the radiative diffusion time at breakout is less than 
the light travel time across the star. The largest red supergiants have a radius of ~ 10^^ 



cm 



(Levesque et al. 2005) so that the longest time of the breakout event is about 1 hour. The 



short timescale of the bursts makes it difficult to detect them, other than with a wide-field 
X-ray telescope. The breakout detection of SN 2008D with Swift ( Soderberg et al.||2008 ) was 
fortunate. 

The situation changes if there is dense mass loss prior to the supernova that creates 
an optically thick region. Calculations of this case go back to early computer simulations 
of supernova light curves. Model 5 of Grassberg et al. (1971) and Model B of Falk & 



Arnett (1977) include a dense circumstellar shell with radius ~ 10^^ cm that determines 
the properties of the shock breakout. The peak luminosity occurs on a timescale of ~ 15 



days in the model of Grassberg et al. (1971). Grasberg & Nadyozhin (1987) considered a 



steady wind of limited duration just before the supernova. Chugai et al. (2004) calculated 



a model for SN 1994W which allows for dense mass loss around the supernova, leading to a 
peak luminosity at an age of ~ 20 days. A systematic numerical study of such events has 



recently been carried out by Moriya et al. (2010). Because of the high surrounding density, 



the supernova can be especially bright and this type of model has been invoked for luminous 
supernovae. In addition to SN 1994W (Chugai et al. 2004), application of the dense mass 



loss model has been made to SN 2006gy (Smith & McCray 2007), PTF 09uj (Ofek et al. 



2010), and SN 2009kf (Moriya et al. 2010), among others 



The aim here is to give an analytical description of the shock breakout process for the 
case where a radiation dominated shock propagates into the mass loss region. Progress on 



this front has been made by Ofek et al. (2010). The model is presented in § 2 and compared 
to observations in 5 3. 



2. INTERACTION WITH DENSE MASS LOSS 



We assume that the dense mass loss can be described by a steady wind acting for 
some time tmi before the supernova explosion; the extent is = v^tmi, where is the 
wind velocity. The actual case may be more complex, but is not well determined; the wind 



assumption has been made in other work (Grasberg & Nadyozhin 1987 Ofek et al. 2010 



Moriya et al. 2010). If the mass loss is in a steady wind, the density = M /Anr'^Vw = 
Dr~^ can be specified by a density parameter, D^:, scaled to a M = 10~^ Mq yr~^ and 
= 10 km s~^ wind so that = 5.0 x lO^^D^r"^ in cgs units. For typical supernova 
parameters, the explosion drives a radiation dominated shock through the star. A radiation 
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dominated shock in a uniform medium has a characteristic thickness described by an optical 
depth Tsh ~ c/vsh, where Vsh is the shock velocity (Weaver 1976). If the wind optical depth 
Tw < c/vsh, then a radiation dominated shock does not form in the wind and the radiation 



diffuses out through an optically thick wind (e.g., Nakar & Sari 2010). In the case that 
Tw > c/vsh, the radiation dominated shock front propagates into the mass loss region. This 
is the case considered here. 

After the shock wave passes through the star and the supernova ejecta tend toward free 
expansion, a shocked layer develops in the wind that is bounded by forward and reverse 
shock waves at Rfs and Rrs, respectively. The shock waves are initially energy conserving. 
Once the forward shock front reaches a place where the diffusion time equals the expansion 
time or, equivalently, ~ c/vgh, radiation diffusion becomes important. We have = 
{Rj} — R~})kD, so the diffusion condition can be written as ~ R~t} — R7,}, where 
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Rd = nDvsh/c = 5.7 X lO^^fsM^-D* cm, where VshA, is the shock velocity in units of 10'* km s 
and k is the opacity k, in units of 0.34 cm^ g~*. Two cases are of interest: R^ > Rd so that 
Rfs ~ Rd when diffusion becomes important, and R^ < Rd so that Rfs ~ Rw at that time. 

In the wind, for diffusion over a distance AR, the diffusion time is td = AR'^/Xc = 
KPwAR? / c. Taking AR R and noting = Dr~^, we have td = k.D/c, independent 
of radius (see also Ofek et al. 2010). If the extent of the optically thick region (to Rph = 



min(i?^, kD)) is not much beyond Rd, then kD/c provides an estimate of the diffusion time 
for the radiation to escape. The time is longer if Rph > Rd- If the diffusion is assumed to 
proceed in steps with size oc R, the timescale is increased by a term that is logarithmic in 
R. For the approximate results presented here, we neglect the logarithmic dependence and 
take 

td = kD/c = 6.6kD^: day. (1) 

We now examine a specific model for expansion into the dense wind. After the shock 
wave has expanded beyond the stellar radius, the ejecta tend toward free expansion {v = r/t) 
in which there is a steep power law section on the outside and a fiat profile in the inner region 
(e.g., Matzner fc McKee||l999 ). Here we make the approximation that pin oc at velocities 
V < Vt and psn oc at v > Vt- The outer profile is less steep than the profile inferred for the 
asymptotic profile expected at high velocities; the wind interaction causes a flatter part of the 
density distribution to be relevant. We then have Vt = {2E / MeY^"^ and psn = Bt^^(r /t)^"^ , 
where B = AE'^ /SrcMe ( Chevalier fc Fransson||l994 ). The interaction between this power law 
density distribution and a surrounding wind, with density p^ = Dr~^, can be described by 
a self-similar solution with adiabatic index 7 = 4/3 (Chevalier 1983), which is appropriate 



for a radiation dominated gas. The radius of the contact discontinuity is 



R 



■cd 



(0.235/D)0-2tO-s 



5.6 X IO^^EsVM-iO-^D; 



o-^to-s cm. 



(2) 
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where £"51 is the explosion energy in units of 10^^ ergs, M^i is the ejecta mass in units of 
10 Mq, and ti is the age in units of 10 days. The highest velocities occur just inside the reverse 
shock wave, where the velocity is Vrs = 0.978Rcd/t = 6.3 x 10^E^^^M-°-^D-^-\°-^ km s-\ 
The applicability of the solution requires that vt < Vrs, or ElfM-^'^D^ti < 32. 

We initially assume that Rd < Rw At time td, 

R,d = Rd = 4.0 X 101^A;°-«E5VM-i°-^D°-^ cm. (3) 
The amount of energy in the interaction shell above some velocity Vq is 

Esh = / -psnV A-Kr dr = ^. (4) 



The internal energy, which is in the form of radiation, is Ej-ad = 0-<i2Esh (Chevalier 1983). 
Using Wo = Rcd(td)/td at time td, we have 

Erad = 0.44 X 10^°A;°-^E5YM-i°-^D°-s ergs (5) 

for the radiative energy in the breakout shell. The timescale for the radiation to escape is 
td, so the luminosity at shock breakout is 

Lb ~ Erad/td = 7.6 X 10^'k-'-^El,^M;,'-'D;^-' erg s-^ (6) 

Taking the radiation energy E^ad and the shocked volume of Atc^R^j-^ — -Rrs)/3 at time td, 
where Rfs and Rrs are the forward and reverse shock radii, the mean radiation energy density 
gives a radiation temperature 

Trad = 6.4 X 10^fc-°-25t-0-25 = 7.1 x 10^ A;-"-^/};^-'^ K, (7) 



where we have used Rfs/Rcd = 1.208 and Rrs/Rcd = 0.978 (Chevalier 1983), and tdi is 
td in units of 10 days. Another estimate of Trad can be obtained from the forward shock 
condition on the pressure: aT^/3 = (6/7)p^f|^ for a 7 = 4/3 gas, where a is the radiation 



constant (see also Ofek et al. 2010). In this case, the temperature coefficient in equation ([T]) 



is 7% higher; the difference is due to the fact that there is a drop in the pressure behind 
the forward shock front in the shocked region. The results assume a smooth wind with an 
density profile, but the temperature of the breakout shell is not strongly dependent on 
this assumption. More generally, if AR is the size of the breakout region, the optical depth 
through the region is pnAR ^ c/vsh and the shock crosses the shell in the diffusion time, 
or AR ^ Vshtd- The resulting density p = c/nvlj^td, when combined with the shock jump 
condition for the radiation pressure, gives a temperature comparable to that in equation ([T]). 
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These results can be compared to the numerical results of Moriya et al. (2010). Taking 
their model sl3w2r20m2e3 with E^i = 3, M^i = 1.3, and D^, = 1, the breakout radius 
in our model is Rd = 5.9 x 10^^ cm, which is inside the outer radius of 2 x 10^^ cm in 
the numerical model. Substituting the parameters into equation (5), the radiated energy is 
Erad = 1-4 X 10^*^ ergs, which is close to the 2.0 x 10^*^ ergs found in the numerical model 



(Moriya et al. 2010). The numerical result may be larger, in part, because the shock wave 



generates power in the more extended circumstellar medium. We also examined the scaling 
of Erad with the parameters and found reasonable agreement with the numerical results. 
The scaling depends on the supernova density gradient, which is only approximately treated 
here. 

Estimates of the observable color temperature of the radiation require considerations 



of whether radiation equilibrium is attained in the emitting region. Following Nakar fc Sari 
(2010), we define a thermal coupling coefficient r] = nBB/(tdhphjf{TBB)), where ubb ~ 



aT^g/3kBTBB is the photon number density in thermal equilibrium, ks is Boltzmann's 
constant, and nph fjiTsB) = 3.5 x lO^^p^T"^/^ s~^ cm~^ is the production rate of photons by 
the free-free process. If sufficient photons are produced to maintain the blackbody number 
density, or 77 < 1, thermal equilibrium is achieved. Using equations ([T| and ([T]) for td and 
Tbb, and p = 7pw{Rfs) (taking into account the factor 7 compression in the shock wave), 
we estimate r] for the breakout shell: 



-0.8 : 



^1.48 



For the standard parameters, the breakout shell is marginally in thermal equilibrium. As 
the radiation propagates into the unshocked mass loss region, the lower density results in 
a deviation from thermal equilibrium. The frequency dependence of the opacity can play a 



role (Moriya et al. 2010) and we do not treat the details of spectrum production here. 



The loss of radiative energy from the shocked region results in the formation of a dense 



shell at radius i?, as seen in numerical simulations (Grassberg et al. 1971 Falk & Arnett 



1977). The expansion of the shell into additional mass loss produces continuing power for the 



supernova, L 
radiative acceleration 



where the wind velocity v^, may be affected by preshock 



The simulations of Moriya et al. (2010) show some evidence for 
acceleration, but it is only significant near the breakout radius because of the dependence 
of the radiative flux, and we neglect it here. The expansion of R can be described by the 
thin shell approximation (Chevalier 1982), yielding R = 0.94Rcd- The resulting power is 



L = 7.1 X 10^3^5\2^-i°-6D°-^tr°-' erg s'^ 



(9) 



The magnitude of the luminosity is similar to that produced by the initial breakout radiation, 
as is seen in numerical simulations (Grasberg & Nadyozhin 1987 Moriya et al. 2010). The 
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luminosity lasts until the shock wave at R reaches the edge of the dense wind, i?^, at 



= l.lEg-f ^M°^252^o.25^L25 ^YieiB R^iq IS in units of 10^^ cm. Fig. la illustrates the 
luminosity evolution with the late flattening from the shell interaction. 



In the case of shock breakout from a red supergiant, the shock front takes ~ 1 day to 
traverse the star and the time for shock breakout is ~ 10^ s (Klein & Chevalier 1978). The 



shock breakout timescale is much less than the time since explosion. In the dense mass loss 
case considered here, the time for the shock front to move to the breakout region is ~ Rd/vgh, 
which is also the timescale for the breakout event. This property of the luminosity evolution 
can be seen in simulations of such events ( Grassberg et al.||1971 Falk &: Arnett||1977 Chugai 



et al. 2004 Moriya et al. 2010). The rise to maximum light can have complications due 



to variations in the gas opacity. At the high circumstellar densities considered here, it is 
likely that the gas is initially neutral and that most of the opacity is due to dust in the 
presupernova environment. The radiation dominated shock wave from the supernova has a 
precursor in the mass loss region that is expected to heat the circumstellar dust. As the 
temperature rises through 1000 — 2000 K, the dust evaporates, giving a decrease in the 
opacity and the photospheric radius drops to where there is a sharp gradient in the opacity 
as the gas becomes ionized. In Type IIP supernovae, this property of the opacity causes 
a recombination wave to back into the expanding envelope with a constant photospheric 
temperature T ~ 5000 — 6000 K. Here, the process is inverted and the photosphere is 
expected to follow the ionization wave moving out through the dense circumstellar gas. This 



phase of approximately constant temperature can be seen in simulations (Grassberg et al. 



1971). Once the circumstellar mass is ionized, the photospheric expansion slows and the 



temperature rises. The light curve rises fairly sharply due to the temperature rise until the 
maximum temperature, and luminosity, is reached. 

We now consider the case that R^ < Rd- The shock breakout process begins at tb ~ 
Rw/vsh- The rise time for the light curve is tr ~ 6R/vsh, where 6R is the distance in 
from Rw across which the diffusion time equals the shock crossing time. We thus have 
{SR)^K,DR~'^/c = 6R/vsh, leading to tr ~ {Rw/vsh){,Rwl Rd)- In this case the rise time can 
be considerably shorter that the time for initial heating of the envelope (Fig. lb). Since 
SR/Ruj < 1, the region involved in shock breakout is not radially extended, so the pressure of 



the escaping radiation can accelerate the gas out to Rw (Ensman 1994). High gas velocities 



are expected around the time of maximum luminosity. A dense shell forms when the radiation 
can escape, but does not produce continuing high luminosity because the dense gas does not 
extend far beyond the breakout point and it has been radiatively accelerated. 



For our specific supernova model, the shock breakout begins when the forward shock 
{Rfs) reaches R^, which occurs at tyj = 16Rl;,ll,E^^'^M^f^D'^-'^^ days, where R^w is R^ 
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in units of 10^^ cm. 



103i?-?;,25^0.5^-0.25^ 



^wlb -^51 ^^-'el 
-0.5£)0.5 



ergs. 



The free expansion velocity at the reverse shock is Vq 
^ 0.25 which leads to a radiated energy of Erad 

The rise time for the light curve is tr = t^Rw/Rd = 41/c~° 

/j. 1 o , , 1 n4,S;_n.S t-)-1.75 

rad 



5.7 X 



0.9 



D2.25 rn 



^a45 25-0.35 ^^yg^ gQ ^j^g^t ^ ^ Erad/tr = 1-8 X lO^^^o-^i^-iY^EifM-i^-^^DO-ss erg s-\ The 
temperature in the shocked shell is like that given in the first part of equation ([T]) except 
that td is replaced by t^,. The temperature is lowered in the escape out to the photosphere 
if there is sufficient photon production in this region (Nakar & Sari 2010). 



3. COMPARISON WITH OBSERVATIONS 



The radiated energy from SN 2006gy from the initial optical rise to around the peak is 



1 X 10^^ ergs (Ofek et al. 2007 Smith et al. 2007), which makes it a good candidate for 



the physical situation described here (with wind optical depth > c/vgh)- In this case, the 
observed rise time gives an estimate of td, the diffusion time. The observations indicate a rise 
time of 60 days. Using equation ([T]), the indicated wind density is ~ 10. The observed 
peak luminosity has sensitivity to the supernova energy. At maximum light the observed 
luminosity was 4 x 10^^ erg s~^ (Smith et al. 2010), which implies E^i ^ 3 for the other 



standard parameters. With these parameters, equation (^ gives Rd = 2.5 x lO^^M^j^"'^ cm 



while the dense medium may extend to ~ 1 x 10 cm (Smith et al. 2010), so that R^ > Rd- 
We attribute the flattening of the observed light curve (Smith et al. 2010) to the continuing 
interaction in the extended region. With i?^ = 10^^ cm the implied circumstellar mass is 



30 Mq. These parameters are close to those found by Smith & McCray (2007) and Smith 



et al. (2010), which is expected because the basic physical picture is similar in the two cases. 
Depending on the ejecta mass, the condition that Vt < Vrs may be violated, but we do not 
expect the results to be strongly affected. 



Smith et al. (2007, 2010) estimated an explosion date of 20 August, 2006 for SN 2006gy; 



max to LjYiaxi 



this time is just before the beginning of a rise in optical luminosity from O.OIL 
where Lmax is the maximum optical luminosity. In the shock breakout view, the sharp rise 
of optical radiation occurs after a time Rd/vsh during which the shock is traveling in the 
optically thick region (Fig. 1). The explosion date is thus ~ 60 days earlier than the estimate 



of Smith et al. (2007, 2010), and mean velocities of uniformly expanding ejecta are lower 



than the estimates of Smith et al. 

Another aspect of the rise to maximum in shock breakout is that the increasing tem- 
perature as the shock breaks out is an important component of the rising luminosity. In the 
case of SN 2006gy, [Smith et al. ( [2010 ) find a temperature T !^ 11, 000 K on days 65 and 71 
(from Aug 20, 2006) (their Fig. 4), which is close to the time of optical maximum light. On 
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day 36, [Smith et al.| (|20To[ estimate T ^ 9, 500 K, usiner the same extinction value as that 
used for other epochs (their Fig. 4). However, they find that with an assumed larger extinc- 
tion, a photosphere with T ^ 15, 000 K provides a better fit to the observed spectrum. The 
higher temperature would bring the temperature evolution more in line with that observed 



in other supernovae, i.e. a decreasing temperature with time, and is advocated by Smith et 



al. (2010). However, in the shock breakout view, the increasing temperature is expected and 



there is no need for a time dependent extinction. 

In our model, the dense shell and shock wave are deep within the circumstellar envelope, 
outside of which is the last equilibrium shell where the spectrum is formed. The radius of this 
shell is not the blackbody radius Rbb because of the nonequilibrium conditions in the medium. 
Outside of the shell is an electron scattering region of moderate optical depth (between c/vgh 
and 1) where the peaked Ha profile with a broad base can be formed ( |Chugai||2001| [Smith et 



al. 2010); the broadening is due to scattering in the thermal gas as opposed to the Doppler 



effect. 



An example of a different type of luminous supernova is SN 2010gx and related objects 



(Pastorello et al. 2010 Quimby et al. 2009). In this case, there is a sharp drop from peak 



luminosity without a flattening, broad lines of O II are present at maximum light, and the 
narrow H or He lines often observed in Type Iln supernovae are not present. The evidence 
points to an explosion in a more compact circumstellar medium for this class of objects. The 
rise time, ~ 50 days, and peak luminosity, (3 — 4) x 10"^^ erg s~^, are comparable to the case 
of SN 2006gy, while the temperature is higher, 15, 000 K. The rise time suggests that is 
not much less than Rd, so R^ ~ Rd- The density is comparable to SN 2006gy. The higher 
temperature in this case can be attributed to the lower opacity due to the lack of H and He, 
and the smaller extent of the opaque circumstellar medium. 

Although the case for dense mass loss years before the supernova appears good, the 



cause of the mass loss is not known. Woosley et al. (2007) suggested that SN 2006gy was the 



result of pair instability eruption before the supernova. Since the radiated energy produced 
in this case is ~ 10^° ergs and the radiated energy in observed luminous supernovae is > 10^^ 
ergs, we have not specifically treated that case here, although if the energy were larger similar 
physical arguments would presumably apply. The dense mass loss is typically attributed to 



LBV (luminous blue variable) eruptions (Smith & McCray 2007), although such eruptions 
are not understood in the context of stellar evolution. An eruption was observed 2 years 



before the explosion of SN 2006jc, which was a H-poor, Type lb supernova (Pastorello et al. 



20071 Foley etaL][2007l ) . The mass loss in this case was too weak to produce the high density 
phenomena discussed here, but it does show the possibility of a presupernova outburst, even 
in the H-poor case. 
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Note added: A model similar to that developed here was presented by |Balberg &: Loeb 



(2011) and applied to the Type lb SN 2008D. 



We thank Claes Fransson for discussions. This research was supported in part by NSF 
grant AST-0807727. 



REFERENCES 



Balberg, S., & Loeb, A. 2011, MNRAS, submitted ( |arXiv:1101.1489[ ) 

Chevalier, R. A. 1982, ApJ, 259, 302 

Chevalier, R. A. 1983, ApJ, 272, 765 

Chevalier, R. A., & Fransson, C. 1994, ApJ, 420, 268 

Chugai, N. N. 2001, MNRAS, 326, 1448 

Chugai, N. N., et al. 2004, MNRAS, 352, 1213 

Ensman, L. 1994, ApJ, 424, 275 

Ensman, L., & Burrows, A. 1992, ApJ, 393, 742 

Falk, S. W. 1978, ApJ, 225, L133 

Falk, S. W., & Arnett, W. D. 1977, ApJS, 33, 515 

Foley, R. J., Smith, N., Ganeshalingam, M., Li, W., Chornock, R., & Filippenko, A. V. 2007, 
ApJ, 657, L105 

Grassberg, E. K., Imshennik, V. S., & Nadyozhin, D. K. 1971, Ap&SS, 10, 28 
Grasberg, E. K., & Nadyozhin, D. K. 1987, AZh, 64, 1199 
Klein, R. I., & Chevalier, R. A. 1978, ApJ, 223, L109 

Levesque, E. M., Massey, P., Olsen, K. A. G., Plez, B., Josselin, E., Maeder, A., & Meynet, 
G. 2005, ApJ, 628, 973 

Matzner, C. D., & McKee, C. F. 1999, ApJ, 510, 379 



Moriya, T., Tominaga, N., Blinnikov, S. I., Baklanov, P. V., k Sorokina, E. I. 2010, MNRAS, 
submitted ( |arXiv :1009.5799|) 

Nakar, E., & Sari, R. 2010, ApJ, 725, 904 

Ofek, E. O., et al. 2007, ApJ, 659, L13 

Ofek, E. O., et al. 2010, ApJ, 724, 1396 

Pastorello, A., et al. 2007, Nature, 447, 829 

Pastorello, A., et al. 2010, ApJ, 724, L16 



Quimby, R. M., et al. 2009, Nature, submitted ( |arXiv:0910.0059| ) 
Smith, N., & McCray, R. 2007, ApJ, 671, L17 
Smith, N., et al. 2007, ApJ, 666, 1116 

Smith, N., Chornock, R., Silverman, J. M., Filippenko, A. V., & Foley, R. J. 2010, ApJ, 
Soderberg, A. M., et al. 2008, Nature, 453, 469 
Weaver, T. A. 1976, ApJS, 32, 233 

Woosley, S. E., Blinnikov, S., & Heger, A. 2007, Nature, 450, 390 



This preprint was prepared with the AAS I^TjrpC macros v5.2. 



- 11 - 




Fig. 1. — Luminosity curves of two types of supernova interaction with dense mass loss: a) 
wind extent greater than the characteristic diffusion radius and b) R^j < R^. In each 
case, there is a time from the explosion at t = to the shock wave reaching a place where 
the radiation can escape and the luminosity rises. In the case R^ > Rd, there is a later, 
slower luminosity decline due to continued interaction of the shock wave (velocity Vsh) with 
slow wind material. 



